Back-to-back correlations of high hadrons in relativistic heavy ion collisions 
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We investigate the suppression factor and the azimuthal correlation function for high pt hadrons in 
central Au+Au collisions at t/s~nn = 200 GeV by using a dynamical model in which hydrodynamics 
is combined with explicitly traveling jets. We study the effects of parton energy loss in a hot medium, 
intrinsic &t of partons in a nucleus, and p± broadening of jets on the back-to-back correlations of 
high pt hadrons. Parton energy loss is found to be a dominant effect on the reduction of the 
away-side peaks in the correlation function. 

PACS numbers: 24.85.+p,25.75.-q, 24.10.Nz 



Recently, from the two-particle azimuthal correlation 
measurements in central Au+Au collisions at the Rela- 
tivistic Heavy Ion Collider (RHIC), the magnitude of the 
away-side jets is found to be significantly suppressed in 
comparison with the near-side jets at high px [lj- In ad- 
dition to the suppression of theyields for high px hadrons 
already discovered at RHIC @, , which was predicted 
some time ago , the high j>x correlation measurement 
may provide a novel opportunity to study the dense QCD 
matter produced in high energy nuclear collisions. This 
phenomenon was first predicted by Bjorkcn in many years 
ago p|- Later, the acoplanarity (transverse momentum 
imbalance) of jets was discussed as a possible diagnostic 
tool for studying the quark gluon plasma (QGP) [f|. 

In this Letter, we investigate the effects of (a) parton 
energy loss during propagation through an expanding hot 
medium, (b) intrinsic fcx of partons in a nucleus, and 
(c) p± broadening of jets on back-to-back correlations 
in Au+Au collisions based on incoherent production of 
mini-jets by employing the hydro+jct model Q. 

The hydro+jet model is mainly composed of two mod- 
els. One is a full three dimensional hydrodynamic model 
which describes the space-time evolution of thermalizcd 
partonic/hadronic matter The other describes the 
dynamics (production, propagation, and fragmentation) 
of high px partons which interact with thermalizcd par- 
tonic matter through a phenomenological model for par- 
ton energy loss. A novel hydrodynamic calculation with 
the early chemical freeze-out shows that the px slope for 
pions is insensitive to thermal freeze-out temperature T th 
P and that the pt spectrum starts to deviate from the 
data at px = 1-5-2.0 GcV/c [ljj. This result naturally 
leads us to include mini-jet components into hydrody- 
namic results. Thus we simulate the space-time evolution 
of both a fluid and mini-jets simultaneously. 

Initial parameters in the hydrodynamic model are so 
chosen as to reproduce dN^/drj in Au+Au collisions at 
^/JW^ = 200 GeV observed by the BRAHMS Collabo- 
ration The transverse profile of the initial energy 
density is assumed to scale with the number of binary 
collisions F° r 5% central collisions, we choose the 
impact parameter as b = 2.0 fm and the maximum ini- 



tial energy density at the initial time tq = 0.6 fm/c as 
39.0 GeV/fm 3 . For further details on initialization in our 
hydrodynamic model, see Ref. 0. We note that contri- 
bution from hard components (defined by the particles 
that have a transverse momentum greater than 2 GeV / c 
just after initial collisions) becomes the order of unity 
and can be neglected when we tune initial parameters in 
the hydrodynamic model. 

We include hard partons using the perturbative QCD 
(pQCD) parton model. The momentum distribution of 
hard partons per hard collision before traversing hot mat- 
ter may be written as 

E ^ p ±=K J2 J g(kT, a )dk 2 T ^g(k T ^)dkl h 



a . b 



da„ 



x J f a (xi,Q 2 )dx 1 f b (x2,Q 2 )dx 2 E^, (1) 

where / and g are the collinear and transverse parts of 
parton distribution functions, x is the fraction of the 
longitudinal momentum and kj> is the intrinsic trans- 
verse momentum of initial partons inside hadron. The 
parton distribution functions f a (x,Q 2 ) are taken to be 
CTEQ5L [13 with the scale choice Q 2 = (px/2) 2 for the 
evaluation of parton distribution. The minimum momen- 
tum transfer pr.min — 2.0 GeV/c is assumed. The sum- 
mation runs over all parton species and relevant leading 
order QCD processes are included. Gaussian intrinsic fcx 
distribution inside hadron with the width of (fc T ) = 1 
GeV 2 /c 2 is assigned to the shower initiator in the QCD 
hard 2 — ► 2 processes. In the actual calculation, we use 
PYTHIA 6.2 [l4| to simulate each hard scattering and the 
initial and final states radiation. In order to convert hard 
partons into hadrons, we use an independent fragmenta- 
tion model in PYTHIA after hydrodynamic simulations. 
A factor K is used for the higher order corrections. We 
have checked that this hadronization model with K = 2.5 
provides good agreement with the transverse spectra of 
charged hadrons in pp collisions |l5j | and of neutral pions 
in pp collisions 01 above pr = 1 GeV/c at y 7 ? = 200 
GeV. The number of hard partons is assumed to scale 
with the number of hard scattering which is estimated 
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FIG. 1: Comparison of Raa(pt) for neutral pions in Au+Au 
10% central collision at y / sjvjv = 200 GeV. A parameter C in 
Eq. @ is chosen as 0.1 (dotted), 0.27 (thick solid), and 1.0 
(dashed) and {k^} = 1 GeV 2 /c 2 is used. We also show the 
hydr ody namic component (thin solid). Data is taken from 
Ref. H|. 

by using Woods-Saxon nuclear density. This assump- 
tion is consistent with the peripheral Au+Au collision at 
y/JW^ = 130 GeV 0. The initial transverse coordinate 
of a parton is specified by the number of binary colli- 
sion distribution for two Woods-Saxon distributions. We 
neglect the nuclear shadowing effect because its effect is 
small at midrapidity at RHIC energies [l7i ITfij. 

Over the past years, a lot of work has been devoted 
to study the propagation of jets through QCD matter 
[ill l2fl l2ll I22Y Here we employ the first order for- 
mula in opacity expansion from the reaction operator 
approach )2'A |23j . The opacity expansion is relevant for 
the realistic heavy ion reactions where the number of jet 
scatterings is small. The energy loss formula for coher- 
ent scatterings in matter has been applied to analyses 
of heavy ion reactions taking into account the expansion 
of the system (2^, |^. The first order formula in this 
approach is written as 

AE = C jf° d T P (T,x(T))(T-T )]n(j^) . (2) 

Here C is an adjustable parameter |2{| and p(r, r) is a 
thermalized parton density, x (t) and E are the position 
and the initial energy of a jet, respectively. In the hy- 
dro+jet approach, p is obtained by using hydrodynamic 
simulations. We take a typical screening scale /x = 0.5 
GeV |2^| and effective path length L = 3 fm. Feedback 
of the energy to fluid elements in central collisions was 
found to be about 2% of the total fluid energy. Hence we 
can safely neglect its effect on hydrodynamic evolution 
in the case of the appropriate amount of energy loss. 

The jet quenching is quantified by the ratio of the par- 
ticle yield in A + A collisions to the one in p + p colli- 
sions scaled up by the number of binary collisions Raa = 



FIG. 2: Energy loss strength dependence on C2(A</>) for 
charged particles in central 0-5% Au+Au collision at ^snn = 
200 GeV. Parameters C in Eq. © are (dotted), 0.1 
(dashed), 0.27 (solid), and 1.0 (dash-dotted). (fc|) = 1 
GeV 2 /c 2 is used. Data is taken from Ref. QJ. 

(d 2 N A+A /dp T dr 1 )/{N hin ^ y d 2 NP+P/dp T dT 1 ). The trans- 
verse momentum dependence of Raa can provide infor- 
mation about the mechanism of jet quenching 0,H3|. In 
Fig. ^ we represent the results of Raa for neutral pions 
in Au+Au 10% central collisions at -^/saw = 200 GeV 
from the hydro+jet model for parameters of C = 0.1, 
0.27, and 1.0. Here we choose b = 2.8 fm for centrality of 
10%. The suppression factor for the model with C = 0.27 
quantitatively agrees with the PHENIX data. 

We now study the dynamical effects of a hot medium 
on the strength of away-side peaks in the azimuthal 
correlation function. The azimuthal correlation func- 
tions of high pt charged hadrons in midrapidity region 
(| rj |< 0.7) are measured by the STAR Collaboration 0. 
Charged hadrons in 4 < pt, trigger < 6 GcV/c and in 
2 < pr.associate < pr,triggcr GeV/c are defined to be trig- 
ger particles and associated particles respectively. The 
relative azimuthal angle distribution (per trigger parti- 
cle) of high px hadrons is written as 

1 f 1A dN 
C 2 (Ac/>) = - / dA v — — — . (3) 

Atriggcr 7-1.4 dAcpdAr) 

Here A(f> and Arj are, respectively, the relative azimuthal 
angle and pseudorapidity between a trigger particle and 
an associated particle. It is found that, in central Au+Au 
collisions (0-5%), the near-side peak (| A(j> \< 0.75 radi- 
ans) has the same strength as the one in pp collisions at 
= 200 GeV and the away-side peaks (| A(j> |> 2.24 
radians) disappear clearly 0]. 

In our theoretical calculations, we take into account 
only contribution from associated particles which origi- 
nate from the same hard scattering as a trigger particle. 
Thus, we compare our results with the STAR data in 
which elliptic flow components are subtracted assuming 
B[l + 2v\ cos(2A^)] with B = 1.442 and v 2 = 0.07 Q. 
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FIG. 3: Intrinsic &t dependence of azimuthal correlation 
function C2(A<j!>). Parameter sets are ((fc 2 ),C) = (1.0,0.27) 
(solid), (2.0,0.35) (dashed), and (4.0,0.5) (dotted). Here the 
intrinsic transverse momenta (k^) are in the unit of GeV 2 /c 2 
and the corresponding energy loss parameters C are chosen 
to fit Raa- Data is taken from Ref. QJ. 



FIG. 4: The effect of p± broadening on azimuthal correlation 
functions C2(A0). Solid line represents the result for C = 
0.35 and (A;*) = 2.0 GeV 2 /c 2 without p± broadening. Dashed 
(dotted) line corresponds to the broadening of model 1 (2). 
See text for details. Data is taken from Ref. Q]. 



As demonstrated in Fig. [3 the away-side peaks are 
gradually decreasing with increasing the amount of par- 
ton energy loss. Nevertheless, we cannot obtain the dis- 
appearance of the away-side peak by the model with 
C = 0.27 which quantitatively reproduces Raa (pt ) a s 
shown in Fig. ^ Only when we underestimate the sup- 
pression factor Raa of the PHENIX data by choosing 
C = 1.0 as depicted in Fig. ^ the away-side peaks 
almost disappear. This suggests that another mecha- 
nism is needed to interpret both the suppression factor 
Raa ~ 0.25 and the disappearance of back-to-back cor- 
relation simultaneously. 

Is a large intrinsic k^ responsible for the disappear- 
ance of back-to-back correlation? If colliding partons 
have only longitudinal momenta directed to the beam 
axis, produced partons through 2^2 processes have the 
same transverse momentum due to kinematics. When 
partons in a nuclcon initially have the transverse mo- 
mentum (intrinsic fcx), this releases the above constraint 
and, consequently affects the back-to-back correlation of 
final high px hadrons. Moreover, multiple interactions 
within the colliding nuclei lead to the so-called Cronin 
effect [2^. This effect can be parametrized phenomeno- 
logically by increasin g th e average intrinsic transverse 
momentum [23.1301 I3TL l32l| . We have chosen in 

the previous analyses (k^) = 1 GeV 2 /c 2 which leads to 
reproduce the pp and pp data at y/s = 200 GcV. To see 
the nuclear effect of intrinsic k-j clearly, we calculate the 
correlation function by choosing (k^) = 2 or 4 GeV 2 /c 2 . 
Here we readjust a parameter C in Eq. (J2J to quanti- 
tatively reproduce Raa(pt) similar quality to (k\) = 1 
GeV 2 /c 2 below 5 GeV/c: C = 0.35 (0.5) for (fc 2 ) = 2 (4) 
GeV 2 /c 2 . Only when (fc£) = 4 GcV 2 /c 2 which is consid- 
erably larger than an estimated value of ~ 2-3 GeV 2 /c 2 



[24l I30I I3H 132^ at fixed target energies, the height of 
the away-side peak becomes small (~ 0.02) as shown in 
Fig. 3. Instead, the near-side peak is also reduced and 
deviates from data. 

Measurements in d+Au collisions should provide an 
important information on the intrinsic fcx at the RHIC 
energy region. The initial fcx broadening should depend 
on the transverse position of the hard scattering. We 
have checked this effect by using the model in Ref. |30l | 
which gives average (k\) ~ 2 GeV 2 /c 2 for central Au+Au 
collisions at RHIC. We found that the effect of the trans- 
verse position dependence of the intrinsic &t on the away- 
side peak is small and conclusion here remains the same. 
It is needed to use the impact parameter dependence of 
the intrinsic fcx for study of the ccntrality dependence of 
the back-to-back correlation. One should fix the intrin- 
sic hr by pp and pA(dA) data in the systematic study 

Mm 

Finally, we take into account the p± broadening of jets. 
Both the non-abelian energy loss and the p± broadening 
of the propagating jets are manifestation of the strong 
final state partonic interactions [ill I20L I2I I2I I3I l34j . 
Partons can obtain a transverse momentum orthogonal 
to its direction of motion in traversing a dense medium. 
As a consequence, traversing partons change their di- 
rection at each time step and follow zig-zag paths in 
the hot medium. The p± broadening of these par- 
tons is quantitatively related to the parton energy loss. 
We compare two models for the average broade ning 
of partons: {p\) = {a s N c / A^dE / dx (model 1) I33 
and (pi) = {a s N c /2)- x C f TQ drp (r, x(t)) (model 2) 
where a s = 0.3 is assumed. Although the formula for 
model 1 is obtained for static plasma, it is accurate ex- 
cept the factor of logarithm in the case of 1+1 D Bjorkcn 
expansion. Here the distribution of p± is assumed to be a 
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Gaussian function. It should be noted that incorporation 
of pj_ broadening does not largely affect the suppression 
factor Raa- In Fig- HQ the correlation functions C2(A0) 
with and without p± broadening are shown for C = 0.35 
and (fc T ) = 2 GeV 2 /c 2 . When the model 1 is used, the 
event average broadening is found to be ~ 2.5 

GeV 2 /c 2 and away-side peaks are suppressed and broad- 
ened. On the other hand, {{p\}) ~ 0.78 GcV 2 /c 2 for 
model 2 and the reduction of away-side peaks becomes 
tiny. This is consistent with the recent study based on a 
pQCD parton model |35j . 

In summary, we have computed the correlation func- 
tion Ca(A^) for high px charged hadrons in central 
Au+Au collisions at v / ijvw = 200 GcV within the hy- 
dro+jet model. We have systematically studied the de- 
pendence of the magnitude of parton energy loss on the 
height of away-side peaks in correlation functions. More- 
over, we also investigated the effects of the intrinsic trans- 
verse momentum of initial partons and the p± broadening 
of traversing jets on back-to-back correlations of high px 



hadrons. About 50% of the reduction of away-side peaks 
comes from the parton energy loss, but its effect is in- 
sufficient to vanish the peaks completely. This indicates 
that the suppression of the away-side peaks occurs by 
complex mechanisms. The intrinsic fcx of partons in a 
nucleus also affects the back-to-back correlation. This is, 
however, due to the increase of energy loss parameter C 
in order to reproduce Raa rather than the increase of fcx 
itself. The p± broadening causes 5-20% reduction of the 
original height of away-side peaks. 
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